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Abstract-  We  describe  dual-mode  ultrasound  phased  arrays 
for  noninvasive  image-guided  surgical  applications.  In  partic¬ 
ular,  we  address  the  problem  of  real-time  visualization  of  ther¬ 
mal  lesion  formation  in  the  target  (e.g.,  tumor)  tissue  using  the 
therapeutic  arrays.  Post  beamfroming  filter  bank  image  recon¬ 
struction  with  nonlinear  compounding  is  utilized  to  improve 
the  lesion  contrast  in  the  (typically)  very  low-contrast  ultra¬ 
sound  images.  It  is  shown  that  the  new  image  reconstruction 
algorithm  leads  to  measurable  improvement  in  lesion  contrast 
on  the  order  of  6  - 15  dB.  This  leads  to  significant  improvement 
in  lesion  detectability  and  size  estimation  by  standard  segmen¬ 
tation  techniques  for  speckle  imagery.  Experimental  results 
strongly  suggest  that  2nd  Harmonic  imaging  could  play  an  im¬ 
portant  role  in  the  enhancement  of  real-time  lesion  visualiza¬ 
tion. 

1.  Introduction 

Image  guidance  has  long  been  recognized  as  the  “enabling 
technology”  for  noninvasive  thermal  surgery.  Highly  re¬ 
fined  energy  application  devices  have  been  developed  and 
for  years.  However,  without  reliable  imaging  techniques  for 
visualization  of  thermal  lesions,  noninvasive  thermal  surgery 
has  failed  to  find  widespread  acceptance  in  the  clinic.  Re¬ 
cently,  image  guidance  methods  based  on  well  established 
imaging  modalities  like  MRI[3],  CT[5],  or  ultrasound  [4,  8] 
have  been  proposed.  Other  imaging  modalities  are  also  be¬ 
ing  developed  and  may  become  available  in  the  foreseeable 
future  [6,  7]. 

An  area  unique  to  ultrasound  that  could  revolutionize 
the  field  of  image  guided  surgery  is  the  development  of  a 
new  generation  of  dual-mode  high-power  phased  array  sys¬ 
tems  capable  of  both  imaging  and  therapy  [8,  9].  These 
piezocomposite  transducers  can  produce  focal  intensity  lev¬ 
els  needed  for  ablative  and  coagulative  thermal  surgery  with 
high  precision.  Furthermore,  the  operating  bandwidth  of 
such  transducers  allows  for  imaging  the  treatment  region 
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with  adequate  image  quality  to  delineate  important  land¬ 
marks  within  and  around  the  target  volume.  With  these  ca¬ 
pabilities,  it  is  possible  to  operate  these  arrays  in  a  “self¬ 
registration”  mode  whereby  the  imaging  capabilities  of  the 
array  are  utilized  in  characterization  of  the  tissue  response 
precisely  at  the  expected  lesion  location.  This  is  due  to  the 
fact  that  the  beamforming  is  common  to  both  the  imag¬ 
ing  and  therapy  modes.  The  main  challenge  to  this  ap¬ 
proach  is  the  low-contrast  nature  of  speckle -ridden  ultra¬ 
sound  images.  This  paper  addresses  a  new  post  beamform¬ 
ing  filter  bank  image  reconstruction  algorithm  with  nonlin¬ 
ear  spatially  weighted  compounding  algorithm  to  mitigate 
this  problem.  The  design  of  this  algorithm  is  motivated  by 
the  nonlinear  nature  of  ultrasound  propagation  in  tissue  and 
the  fact  that  thermal  lesions  are  known  to  have  increased 
level  of  harmonic  generation.  Experimental  results  demon¬ 
strate  the  potential  of  the  new  algorithm  in  both  enhancing 
the  lesion  contrast  and  size/shape  analysis. 


2.  Image  Formation  Model 

Figure  1  summarizes  the  image  acquisition  and  image  for¬ 
mation  model.  A  64-element  array  optimized  for  maximum 
energy  delivery  at  1  MHz  operating  frequency  is  used  for  le¬ 
sion  formation  in  sample  tissue.  Lesions  are  formed  by  fo¬ 
cusing  the  array  at  a  point  within  the  target  and  maintaining 
high-power  output  for  time  intervals  on  the  order  of  seconds 
(1-5  seconds  typical).  The  power  is  interrupted  for  short 
intervals  (milliseconds)  to  acquire  image  data  by  transmit¬ 
ting  short  (/js)  pulses  from  all  64-elements  and  receiving 
on  selected  elements  using  a  matrix  switch.  Once  the  im¬ 
age  data  set  is  collected,  RF  beamforming  is  performed  to 
form  standard  echographic  images  of  the  target  region.  It 
is  well  known  that  these  images  have  very  low  contrast  due 
to  the  coherent  nature  of  image  acquisition  which  produces 
speckle.  This  renders  the  lesion  detection  very  difficult  us¬ 
ing  ultrasound. 
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2.2.  Single -Transmit  Imaging 


Fig.  1.  Image  Acquisition  and  Image  Reconstruction  Algo¬ 
rithm 


P(  Xp,  Zp ) 


Fig.  2.  Coordinate  system  used  in  the  synthetic  aperture 
imaging  system. 

2.1.  Synthetic  aperture  imaging 

Images  were  obtained  using  a  full  synthetic  aperture  tech¬ 
nique  [10].  The  image  pixel  at  coordinates  ( xp,zp )  was 
therefore  computed  by  (Fig.  2): 
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where  i  is  the  transmit  element  index,  j  is  the  receive  el¬ 
ement  index,  A,;  is  the  transmit  apodization  weight  at  el¬ 
ement  i,  Bj  is  the  receive  apodization  weight  at  element 
j.  Rip  and  T!JV  are  the  distances  from  the  transmit  and  re¬ 
ceive  elements,  respectively,  from  the  image  pixel,  c  is  the 
speed  of  sound  in  the  medium  being  imaged,  and  is 

the  echo  acquired  when  transmitting  with  element  i  and  re¬ 
ceiving  with  element  j.  Specialized  image  reconstruction 
programs  were  written  to  evaluate  the  intensity  at  (xp,zp) 
in  either  ( r ,  sin  9)  or  Cartesian  coordinates. 


The  synthetic  aperture  imaging  algorithm  described  above 
can  be  used  to  produce  the  highest  quality  conventional  im¬ 
ages  that  can  be  expected  from  a  given  array.  However,  due 
to  the  fact  that  transmit  beams  are  synthesized  by  super¬ 
position  of  single-element  transmit  patterns,  the  nonlinear 
interactions  of  the  real-time  transmit  beams  cannot  be  ac¬ 
counted  for  using  this  imaging  algorithm.  Therefore,  we 
have  modified  our  64-channel  phased  array  driver  to  allow 
for  pulsed  transmission  on  all  64  channels  simultaneously. 
This  allowed  us  to  use  the  full  power  of  the  transmit  beams 
and,  therefore,  observe  their  nonlinear  interactions  with  the 
tissue  media.  The  image  formation  process  due  to  a  single- 
transmit  focused  beam  is  a  modified  version  of  Eq.  (1)  as 
follows: 
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where  Sj(t )  is  the  received  waveform  at  element  j  due  to 
the  transmitted  beam  and  Ro  is  a  fixed  distance  determined 
by  the  focal  depth  of  the  transmit  beam.  All  other  quantities 
in  Eq.  (2)  are  the  same  as  their  counterparts  in  Eq.  (1). 
It  should  be  noted  that  beamforming  is  performed  in  the 
RF  domain  using  true  delays,  i.e.,  no  baseband  conversion 
is  done.  This  retains  all  the  frequency  components  in  the 
beamformed  data  for  postbeamforming  filtering  operations 
described  below. 

2.3.  Post  Beamforming  Filtering 

Since  beamforming  is  performed  in  the  RF  domain,  all  the 
frequency  components  in  the  received  echo  signals  are  re¬ 
tained.  This  is  important  since  the  echo  signals  can  be  ex¬ 
pected  to  contain  a  mix  of  harmonics  (and  possibly  subhar¬ 
monics)  depending  on  the  nonlinearity  of  the  tissue  medium 
being  imaged.  This  is  especially  true  for  the  single-transmit 
imaging  mode  where  the  transmit  power  is  sufficiently  high 
to  produce  harmonic  components  in  nonlinear  tissue  me¬ 
dia.  Postbeamforming  filtering  can  be  used  to  isolate  spe¬ 
cific  harmonic  components  and/or  enhance  axial  resolution 
if  the  SNR  of  the  system  is  sufficiently  high.  Algorithms  for 
postbeamforming  filterbank  image  reconstruction  for  pulse- 
echo  ultrasound  can  be  found  in  [11].  For  the  purposes  of 
this  paper,  the  filter  bank  is  formed  of  bandpass  filters  with 
center  frequencies  at  a  set  of  preselected  harmonics  (or  sub¬ 
harmonics). 

2.4.  Nonlinear  Compounding 

Ultrasound  propagation  in  tissue  media  is  inherently  nonlin¬ 
ear.  RF  echo  signals  obtained  using  modern  scanners  with 
wideband  transducers  carry  harmonic  components  that  are 


generated  by  the  medium  nonlinearity,  i.e.,  they  are  not  part 
of  the  originally  transmitted  imaging  pulses.  Since  the  tis¬ 
sue  nonlinearity  parameter  exhibits  higher  contrast  between 
various  tissue  compared  to  tissue  reflectivity,  imaging  this 
parameter  can  lead  to  higher  contrast  images  of  soft  tissue. 
For  a  variety  of  physical  reasons,  it  is  well  known  that  ther¬ 
mal  lesions  generate  higher  harmonics  at  levels  higher  than 
normal  tissues.  By  careful  design  of  the  transmit/receive 
beamforming  and  the  post  beamforming  reconstructions  fil¬ 
ters,  high  contrast  images  of  lesions  can  be  formed  from 
ultrasound  images  formed  at  the  fundamental  and  some  of 
its  higher  harmonics.  In  general,  an  optimal  compound¬ 
ing  rule  must  be  derived  based  on  statistical  model  of  the 
imaged  region  and  the  expected  lesion  size/location.  Fur¬ 
thermore,  the  quality  of  the  beamforming  at  different  har¬ 
monics  or  scales  will  be  spatially  heterogeneous.  For  ex¬ 
ample,  since  the  array  sampling  function  (in  wavelengths) 
is  coarser  at  the  higher  harmonics,  images  formed  at  these 
frequency  will  have  high  contrast  along  the  axis  of  the  trans¬ 
mit  beam,  but  the  contrast  deteriorates  away  from  the  main 
axis.  Furthermore,  since  the  harmonics  are  orders  of  mag¬ 
nitudes  smaller  than  the  fundamental,  compounding  is  per¬ 
formed  with  log-compressed  images  after  appropriate  scal¬ 
ing.  Therefore,  with  reference  to  Figure  1,  the  compounding 
is  performed  in  two  stages: 

1.  The  outputs  of  the  filterbank  are  envelope  detected, 
scaled,  and  log  compressed. 

2.  A  spatially-weighted  sum  of  the  harmonic  compo¬ 
nents  is  performed.  The  weighting  function  is  derived 
from  the  spatial  contrast  function  (SCF)  of  the  imag¬ 
ing  array  at  the  different  harmonics. 

Assuming  that  the  envelope-detected,  normalized  and  log- 
compressed  images  obtained  from  the  filterbank  are  given 
by  I\ ,  1-2 ,  •  •  • ,  In,  then  the  compounded  image  is  given  by: 

N 
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where  IT,;  is  a  weighting  function  reflecting  the  relative  en¬ 
ergy  at  the  harmonic  with  respect  to  the  fundamental  at  the 
(expected)  lesion  location  and  Si(x,y)  is  the  SCF  associ¬ 
ated  with  the  /th  harmonic.  The  SCFs  at  the  different  har¬ 
monic  components  are  computed  from  the  inverse  of  the 
beam  patterns  when  the  transmit  beam  is  focused  at  the  cen¬ 
ter  of  the  (expected)  lesion. 

3.  Results 

Images  shown  in  Figure  3  and  Figure  4  demonstrate  the 
basic  ideas  behind  the  proposed  image  compounding  algo¬ 
rithm.  The  figures  show  images  obtained  before  and  after 


lesion  formation  in  an  ex- vivo  liver  tissue  samples.  The  ex¬ 
pected  lesions  are  cigar-shaped  with  length  of  10  mm  and 
width  of  2  -  3  mm.  All  images  are  shown  at  25  dB  dynamic 
range  and  have  the  same  axial  and  lateral  extent.  In  both 
experiments,  the  lesion  was  formed  at  the  geometric  cen¬ 
ter  of  the  array  and  was  expected  between  90  and  100  mm 
in  the  axial  direction.  All  images  show  the  front  surface 
of  the  sample  at  80  mm  and  the  sponge  backing  phantom 
at  120  mm.  For  each  experiment,  two  pairs  of  images  are 
shown,  one  before  and  one  after  the  lesion  formation.  Im¬ 
ages  shown  on  the  left  hand  side  of  each  figure  show  the  tar¬ 
get  before  lesion  formation  while  images  on  the  right  hand 
side  show  the  target  after  lesion  formation.  The  images  on 
top  of  each  figure  are  formed  at  the  fundamental  (1  MHz 
50%  bandwidth)  while  those  at  the  bottom  of  each  figure 
are  formed  at  the  2nd  harmonic  (2MHz  40%  bandwidth). 
Each  pair  of  images  is  displayed  on  a  log  scale  and  normal¬ 
ized  such  that  the  0  dB  corresponds  to  the  maximum  of  the 
image  on  the  RHS,  i.e.,  the  image  of  the  target  after  lesion 
formation.  The  first  experiment  represent  a  fairly  homo- 
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Fig.  3.  Images  before  and  after  formation  of  Lesion  1. 

geneous  target  tissue  while  the  second  has  a  pair  of  blood 
vessels  at  the  target  location  as  can  be  seen  from  the  images 
taken  before  lesion  formation.  While  both  the  fundamen¬ 
tal  and  harmonic  images  show  enhanced  lesion  contrast  in 
both  cases,  the  harmonic  images  show  a  net  increase  in  con¬ 
trast  by  22  dB  in  the  second  case.  On  the  other  hand,  the 
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4.  Conclusions 

We  have  shown  that  post  beamforming  filterbank  recon¬ 
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rently  underway. 
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